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END 



length = 64; //fast fourier transform length ^ 

delta = .1; %input ( ' Enter stop band tolerance: '); 
beta = 0.2; %input ( T Enter pass band tolerance: '); 
wp=8 ; 
ws=7 ; 

mu = 28; %CP length 

load channel_impulse_resp,dat; %load channel impulse response 

N=19; % sirf filter length 

ch = channel__impulse_resp; 

mm=max(abs (ch) ) ; 

ch=2047*ch./mm; 

NN=size (channel_impulse_resp) ; 
M=NN ( 1 ) ; 
for i=l:mu 

mv(i)=l; %initialize mv 

end 

O conv (mv, abs (ch) ) ; 

[Y,i] = max(abs(C)); %find max which is beginning of GI (CP) 



% construct the channel impulse response matrix 
for i=l:M 

for j=l:N 

if i-j+l<=0 

break; 
else 

H(i,j)=ch(i-j+l) ; 

end 

end 

end 
n=l; 

for i=M+l:M+N-l 
n=n+l ; 
k=0; 

for j=n:N 

k=k+l; 

H (i, j ) =ch (M-k+1) ; _ 

FIG. 4A 

end 



for i=l:mu 
\ mv(i)=l; 
\ end 

V C= conv (mv, abs (ch) ) ; 
V [Y, I] = max (abs (C) ) ; 



for i=l:mu 
mv (i) =1; 
end 

C= conv (mv f abs (ch) ) ; 
[Y, I] = max (abs (C) ) ; 



Id 



f % initialize the SIRF filter to an arbitrary value 
c=[ -0.0593 
0. 01047 
-0. 062386 
0.02418 
-0.065114 
0.030031 
-0.039083 
0.13789 
-0.10266 
-0.014681 
-0.11224 
0. 0546 
-0.12642 
0. 18608 
-0.020895 
0.38407 
-0.30117 
-0.37885 
0. 42326] ; 



(~<3 = f ft (c, length) ; 
Sf = abs (g) ; 
7 lamda = 1; 
/ tol = 50; 
V, cold = sum (abs (c) ) ; 
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'%begin iteration 

for ii=l:20 % 20 is the number of iterations 
% projection on the the set C2 
for m = wp: length/2 

if ( f (m) > (1+beta) ) 
gg=g (m) ; 

g(m) = (1 + beta) /f (m) ^complex (real (g (m) ), imag (g (m) )) ; 
g(m) = gg-flamda* (g (m) -gg) ; 
end 

if( f(m) < (1-beta)) 
gg=g(m) ; 

g(m) = (1 - beta) /f (m) ^complex (real (g (m) ), imag (g(m) )) ; 
g(m) = gg+lamda* (g (m) -gg) ; 

end 

end 

for m = l:ws 

if (f (m) > delta) 
gg=g (m) ; 

g (m) = delta/f (m) ^complex (real (g (m) ) , imag (g (m) ) ) ; 
g (m) = gg+lamda* (g (m) -gg) ; 

end 

end 

for m =. 2: length/2 

g(m+length/2) = conj (g (length/2-m-h2 ) ) ; 

end 

cr = real (if ft ( (g) .length) ) ; %transform to time domain 
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time domain projection on C4 set 
for i = 1:N 

c(i) = cr(i); %project into the set C4 

end 

time domain projection on C5 set 
for n=I-mu-10 : I-mu 
norm = 0; 
prod = 0; 
%prodl=0; 
for i=l:N 

norm = norm + H(n,i) ^ 2; 
prod = prod + c(i)*H(n,i); 
%prodl = prodl + c (i) *ch (n-i+1) ; 

end 

if (prod > tol) 
for nn=l:N 

cc = c (nn) ; 

c(nn) = c(nn) + ( (tol-prod) /norm) * H(n,nn); 
c(nn) = cc+lamda* (c (nn) -cc) ; 

end 

end 

if (prod < -tol) 
for nn=l:N 

cc = c(nn) ; 

c(nn) = c(nn) + ( (-tol-prod) /norm) * H(n,nn) ; 
c(nn) = cc+larada^ (c(nn)-cc) ; 

end 

end 
end 

for n=I:I+100 
norm = 0; 
prod = 0; 
for i=l:N 

norm = norm + H(n,i) * 2; 

prod = prod + c(i)*H(n,i); 

end 

if (prod > tol) 
for nn=l:N 

cc = c (nn) ; 

c(nn) = c(nn) + ( (tol-prod) /norm) * H(n,nn); 
c(nn) = cc+lamda* (c (nn) -cc) ; 

end 

end 

if (prod < -tol) 
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^end 



for nn=l:N 

cc = c(nn) ; 

c(nn) = c(nn) + ( (-tol-prod) /norm) * H(n,nn); 
c(nn) = cc+lamda* (c (nn) -cc) ; 

end 

end 

end 

ss = abs (cold-sum(abs (c) ) ) 
cold = sum(abs (c) ) ; 

g = f ft (c, length) ; %transform to frequency domain 
f = abs (g) ; 



for i = 1:N 

sirf(i) - c(i); %sirf will hold the SIRF filter coefficients 
end 
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